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Abstract 

The problem of coordination and control of multiple spacecraft (MS) moving in formation is con- 
sidered. Here, each MS is modeled by a rigid body with fixed center of mass. First, various 
schemes for generating the desired formation patters are discussed, Then, explicit control laws for 
formation-keeping and relative attitude alignment based on nearest neighbor -tracking are derived. 
The necessary data which must be communicated between the MS to achieve effective control are 
examined. The time-domain behavior of the feedback-controlled MS formation for typical low- 
Earth orbits is studied both analytically and via computer simulation. The paper concludes with a 
discussion of the implementation of the derived control laws, and the integration of the MS forma- 
tion coordination and control system with a proposed inter-spacecraft communication/computing 
network. 


l Introduction 

Future use of multiple micro-spacecraft moving in formation for space exploration, constellation 
space antennas and interferometers, and space-based global communication systems calls for novel 
approaches to the design of spacecraft control systems. I n particular in recent years, growing 
emphasis is placed on the concept of separated spacecraft interferometry (SSI). The SSI concept 
envisioned the collecting apertures to be located on separate spacecraft while central combining 
instruments to be located on yet another spacecraft. A virtual structure is therefore developed with- 
bout the real need for maintaining the necessary structural rigidity. The SSI provides measurements 
unachievable with other techniques and allows long baseline lengths and orient ation changes. 

This paper focuses on the development of a control system architecture for the coordination and 
cent rol of a fleet of micro-spacecraft moving in formation. Here, we are dealing with a collection 
of systems which interact with each other in a cooperative manner to achieve a common objective. 
Although control of a single spacecraft is based on well-established control theory concepts and 
methodologies, the control systems for multiple spacecraft moving in formation require architectures 
which differ from those of conventional single spacecraft control systems. To provide a desired 
formation, basic mathematical models for controlled movement of rigid bodies in free space is 
presented. The control laws for the coordination of spacecraft attitude during motion to achieve 
a specified objective (e.g. orient each spacecraft along a given direction) is also developed. This 
is followed by the derivation of control laws for formation keeping and relative attitude alignment, 
The time-domain behavior of the feedback-controlled formation flying for typical low-Earth orbits 
is studied both analytically and via computer simulation. Emphasis is placed on determining the 
information exchange needed for achieving and maintaining a desired formation, and the conditions 
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for ensuring formation stability in the presence of various types of perturbations. Other important 
factors such as collision-avoidance, spacecraft failures, loss of communication between spacecraft, 
and also various constraints imposed by implementation and the physical size of the micro-spacecraft 
are discussed. 


2 Modeling of Multiple Spacecraft in Formation 

To simplify the development, we consider only a spacecraft triad modeled by rigid bodies with 
fixed centers of mass moving in free space under the influence of a gravitational field and external 
disturbances. We introduce the following coordinate systems in the three-dimensional Euclidean 
space 3? 3 : (i) an inertial coordinate system To with orthonormal basis &o = {ex,eY,ez}, and 
(ii) a set of moving coordinate systems T t ,i = 1,2,3 whose origins 0, : are at the mass centers of 
the spacecraft. Let B<, = {ei x ,ei y ,ei Z }denote an orthonormal basis associated with the moving 
coordinate system (abbreviated by “MCS” hereafter) Ti, and [w] 7 the representation of the vector 
w with respect to basis Bj. The the basis vectors in 13 0 and Bi are related by a linear t ransformat ion 
Cj defined by 


&ix — C{y — CjC y , &iz C{C Z , i — 1,2,3 (l) 

whose representation with respect to basis Bq is given by the direction cosine matrix C(qj) = 
(Qii - + ^qaT- 2qi 4 Q(qi), where 


Q{qd = 


0 

-Qi3 

qi 2 

^ A 

; qi = 

Qi 1 

qtz 

0 

-qi i 

qa 

_ -qi.2 

Qi i 

0 


© 3 _ 


( 2 ) 


andq,; = [qf,Qi 4p denotes the unit quaternion with q^ being the Euler symmetric parameters [5] 
defined by qij = SijSin((j)i/2), j = 1,2,3; q l4 = cos{4>i/2), satisfying the constraint = 

1, i= 1,2,3 where (pi is the principal angle and the £;/s are the components of the principal vector 
of rotation £; defined by £, = s t i e x + e-ae y + e l3 e z = e\\e vx + e^e-iy + The time derivitive 

of qij is related to the angular velocity = u ix c ix + u iy eiy + u iz e iz of Ti relative to the inertial 
coordinate system by 


<jq± 

dt 

d q« 


= OJi x qJ/2, 

= -(w», q,)/2, 


(3) 


where w x v and w ■ v denote respectively the cross and scalar products of vectors w and v in 

7Z 3 . 


Let djdt 0 and djdti denote the time derivative operators with respect to and Ti respectively; 

and Vi the time derivative operator d/dt Q in Ti defined by Vi w =dwfdti + u>i x w for w e 73 3 . 
The angular velocities of the i-th spacecraft or T relative to To are given by the following Euler's 
equations relating that time derivative of the angular momentum Ii^i with respect to To to the 
control torque t c ;: 


Vi(IiWi ) = d(IiU>i)/dti + u>i x (IiCOi) = Iidui/dti + w ; x (J 7 ;o>j) = (4) 

for?; = 1,2,3, where Ii is the tensor of inertia associated with the i-th spacecraft. The time 
derivative of w; may also be taken with respect to T 0 . 

In formation acquisition and keeping, we are interested in the relative motion between any pair of 
spacecraft labelled by subscripts i andj. Let the MCS.F' be T as defined earlier (SeeFig. 1), and 
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p,j { denote the position vector of the j-th spacecraft relative to Ei. The evolution of p^ with time 
is governed by 


VfiPji) = ( fcj + fgj)/Mj - (f ci + f gi )/Mi , (5) 

where Vf denotes the time derivative operator d?/dtl in the MCS -T 7 ; given by 

Vf = d?/dt% + (dui/dti) x • -f 2u>i x d/dt; + u>i x (cvj x •). (6) 

3 Formation-Keeping Control Laws 

The movement of spacecraft triad in formation can be achieved in marry ways. Here, we use a 
simple approach proposed in [6], Let the first spacecraft be designated as the leader, and the 
remaining spacecraft as followers. The leader provides the reference motion r\=ri(f) (relative 
to E 0 ) for the followers, and also two nonzero deviation vectors hi (t), i = 2,3, such that the point 
set V(t) ={ri (t), r*i (t) + h 2 (t), r*i(f)+ /13(f)} defines the desired formation pattern at time t. It 
is assumed that the convex hull of P (t) is a simplex. We define the positional error for the i-th 

follower relative to the leader as Ei(t) = hi(t)~ pn(t),i =2,3. 

In what follows, we shall derive formation-keeping control laws assuming that each spacecraft in 
the triad moves in a low Earth orbit (LEO), and the desired motions for the followers are given by 
di(t) =n(t ) +hi(t), i= 2,3. Using (5), we obtain the following differential equation for Ep 

Vf(Ei) — E- t -+- (1){ x Ei -T 2 u>i x Ei + oji x (u x Ei) 

= T^i(hi) + f gi/Mi — f gi /Mi + u c \ — Ud, (7) 

where u C j =f C j/Mj. 

Assuming a central Newtonian gravitational force field, the gravitational force acting on the i-th 
follower has the form: f gi -—pMiri/\\ri\\ 3 , where p is the geocentric gravitational constant; ri 
is the vector specifying the position of the mass center of the i-th spacecraft relative to the inertial 
frame fF 0 \ and ||r,;|j the Euclidean norm of r*. Let e ri =r ! ;/||r,;||. We can write 


r { - ||ri||e r?; 

= ll^i II {(®ri ■ &ix)&ix T (e r j . eiy)eiy + . e,; 2 )ei z }, (8) 

where rj = r,; + Pji ■ For LEO spacecraft triad moving in formation at approximately the same 

altitude, we have ||r*||/||ri|| = I. Thus, the components of {f gl /M\-f gi /Mi) in (7) with respect 
to basis Bi can be approximated by 


(. fgi/Mi - f gl /Mi) ik = <4 0 {\\ri\\{e ri . e ik ) - (||T-ii|/||r<||) 3 (/? U fc + ||n|j(e ri . e ik ))} 

= -UioPlik, k = x,y,z, (9) 

where uif 0 (t)=p/\\ri(t)\\ 3 is the orbital angular speed of the i-th spacecraft about the origin of T 0 at 
the timet. Substituting (9) into (7), making use of the identity u>jx(a>jx Ei) =(^*-£ , j)a;i-||aj.|| 2 J E; j;) 
and assuming u>i 0 = co 0 (a positive constant) lead to 


Ei + oji x Ei 4- 2 ojj, x Ei + (a ■ Ei)u>i 


■ io 
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To derive a formation-keeping control law for the i-th follower, we consider the time rate-of-change 
of the positive definite function V =( KuEi . Ei +Ei-Ei )/2, I< u >0. It can be verified [8] that 
if we set 


u c i ’ I- 1 (-c Oi x + r ci ) x {hi - Ei) + (a >; . ( h, - Ei))u>i 

+ ( W L — ll 2 ) — Ei) + KuEi + K 2 iE{ + 2<jjj x hi + hi + u ci 5 (11) 

where h-i =d,hi/dti,’hi =d?hijdtf, and # 2 ?: is a positive constant, then dVn jdt = —Ku \\E t 11 2 <0. 
The equation for [Ei]i corresponding to the feedback-controlled system is given by 


[Ei]i + {I<2i[I} + 2Q([a> i ] i (t)))[JE? < ] i + KulEili = 0 (12) 

Since Q([wj],;) is slew-symmetric, all solutions [Ei, Ei]i (t ) — ► (O, O)^ 0 as t + oo for any 
Ku,K2i >0 . 

Remarks. Rl : The model -dependent Control law (11) corresponds to state-feedback linearization 
controls which involve partial cancellation of the terms in (10). Assuming perfect cancellation, 
there is no coupling between the equations for the tracking errors of followers given by (12). Since 
perfect cancellation is not achievable physically due to inaccurate knowledge of the model parameter 
values, sensor errors, act uator saturation, and unmodelled external disturbances, it is of import ante 
to determine the effect of imperfect cancellation on the behavior of the feedback-controlled system. 
Here, we model this imperfection by introducing a persistent disturbance N in (12) as follows: 


[Ei} i+ (K 2 i[I} +2Q(M i (t)))[J3 i ] < + KulEili 

= N(t, [wi(t)]i ■ [pji(t)]», [u cj (t)]i, [Ei]i, [i?i]i) s (13) 

and require that the zero state is totally stable [7], Since B\ is a skew-symmetric matrix, the zero 
state of (12) is uniformly asymptotically stable for any Ku, K 21 >0. Then, total stability of the 
zero state follows from a well-known theorem of Malkin [7], 

R2: The control laws (11) for the followers require the knowledge of its own attitude control law r ci 
and the control law u c 1 of the leader. The latter information must be transmitted to the follower 
spacecraft. Note also that control laws (11 ) can be rewritten as 


Uci = Ii 1 {-Ui X (I»Wj) + Td) X Ei - {iOi ■ Ei)u)i - (Ui 0 - ||cj,;|| 2 )jE7,; 

+KuEi + K2i,Ei + U ci + T>f{hi) + u-’i 0 hi. (14) 

The terms u c i and Xf( hi ) +u 2 0 hi in (11') correspond to a feed-forward control. When the norms 
of these terms are large, the norm of «« is also large. This situation may be alleviated by replacing 
the term by a suitable scaling depending on the norm of V\(hi) + u>i 0 hi. I n the important special 
case where the spacecraft move in a nearly circular LEO and the deviation vector hi rotates about 
the Earth's center with angular velocity w* = u,: 0 , then V \ (hi) + Jfjii = 0 or hi is close to a 
solution of the simple harmonic oscillator equation d 2 hi/dtl+uf 0 hi =0. 

R3: It is evident that if each follower applies control law (11), then the desired formation pattern 
V = P(i), t> 0, is asymptotically stable, i.e. given any real number e >0, there exists a 8 > 0 
such that A (t) <6 =$>• A (t) < e for all t >0. Moreover, A(t)— ► O as f->oo, where 
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( 15 ) 


/ 3 \ 1 / 2 

A(i) = EW^II^WII 2 + CT 2*||^(*)II 2 } , 

\i=2 / 

and £7ii,^2t are specified positive weighting coefficients. Here, asymptotic stability is only local in 
the sense that the convergence of A (f) to O as f — > 00 is attained if the deviation of the initial 
formation pattern at t =0 from the desired one is sufficiently small. In physical situations, the 
possibility of collision between spacecraft must also be considered. 

R4: It has been shown recently [9] that under certain mild conditions, asymptotic stability of the 
formation pattern V =V(t) can be achieved using simplified versions of control laws (11') given by 

u d = K\iEj. + I\.2iEi + {u c i + T>f {hi) + cuf 0 hi}, i = 2,3, (16) 

which correspond to proportional -plus-rate feedback plus feed-forward controls. The inclusion of 
the feed-forward control is essential for asymptotic stability. 


4 Attitude Control 

Let the desired attitude and angular velocity of the i-th follower at timet relative to the inertial 
coordinate system T 0 be specified respectively by the MCSJ^(f)and @(t), which may depend 
on the attitude and angular velocity of the leader, eg. the desired attitude and angular velocity 
of the i-th follower correspond exactly to T\(i) and uj i(t) respectively. It is of interest to control 
the relative attitudes and angular velocities between the spacecraft. Here, we shall derive control 
laws for the followers which are expressed in terms of their instantaneous attitudes and angular 
velocities of relative to the inertial coordinate system T 0 or to the MCS^n 

Let the unit quaternion corresponding to Ef(t) relative to the inertial coordinate system T 0 be 

denoted by qf(t) = [g^(t),qf 4 (t)] T . We assume that qf and uf are consistent in the sense that 

Ad 

they satisfy (3) and (4) with" control torque rf. We introduce the deviations 6q { = q\- q L = 
[qf - q ( L - <m] T a n d@ = u>f - It can be verified that Sq t satisfies 

^ = (Qi^i - Quvfuf x cfi +u> ?; x <?;)/ 2. (17) 

= -(“i ■ cfi-Ui, &•)/ 2 . } 

To derive an attitude control law, we consider the following positive definite function Vu =K q iVu + 
Vu defined on R 7 , where I< qi is a given positive constant and Vu = Sqf A + 8q t . 8q { \ Vu = (8u>i . 

Ii8iOi)/2. The time derivative of Vu along the solutions of the equations for 6qi,6ui is given by 

~ - «9i49i - 9? x Sq t ) + 4 - Tc , - x (/,&*)/ 2} ■ Su,,. (18) 

Thus, if we set 

Tci = Kgiiqfyqi - Squqf - qf X 8%) + 7& - u>f X (1,8^^/ 2 + KuiliSuti, (19) 

where K U i is a positive constant, then dVufdt =-K wi 8u)i . 1,8 < 0 and Vii(t) <Vu (0) for all 
t> Oimplying uniform boundedness of ||<5u;{(f)|| for all t>0, By considering d 2 Vu/dt?, and making 
use of Barbalat’s Lemma [10], we can deduce that ( dVu/dt ) ( t)->Oast-> 00, or — > u>f(i) as 

t— mx). B ut it does not follow that 8qn(t)— >0 amlSq^t )— > O ast-+ 00. 
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To proceed further, we make use of the fact that the quaternion (Ag i5 At^) of the desired attitude 
or Ff relative to Ti is related to the quaternions (qf, qf 4 ) for (Ff, and {q^q% 4 ) for T t relative to the 
inertial coordinate system .F 0 by 

= qmt - qi&i -qiXqf, A q iA = q H qf 4 + q { ■ q$. (20) 

When coincides with Tf (i.e. bq,- =0 andSqu = O), we have A q { = O and Aq i4 =1. Using 
(18), control law (17) can be rewritten as 

r c i = -Kgi A <h + r c ;_ uf x (IiSwi )/ 2 + 1^1,8^. (21) 

Following an analysis similar to that given in [ii], conclude that Aq^ and 6ui. t (t ) — ► Oas T — >oo 
for any positive K q i and K u i. 


5 Implementation of Control Laws 

We observe that the implementation of control laws (11) and (19) for formation keeping requires 
a knowledge of [Ei]i,[Ei]i,[^i]i',[r c ili and \ud\j. at any time t. The quantities [E .-] * and [Ej\i can 
be determined from [pii^and [ pn]i which require measurement of the position and velocity of the 
leader relative to the i-th follower spacecraft. These quantities can also be obtained by transmitting 
the position and velocity of leader to the i-th follower. Also, the control of the leader at any time 
must also be transmitted to the i-t h follower spacecraft. When one or more spacecraft failure 
occurs, one may adopt the following backup schemes for control law implement aiton depending on 
the nature of failure: 

(i) Inter-spacecraft Communicaiton System Failure: One may obtain estimates of [pii]i and [pu}i 
by using on-board optical range sensors, or by setting the relative position and velocity between 
the failed and active spacecraft at their nominal values temporarily until the failure is recovered. 

(ii) Overall Spacecraft Failure: Here the failure is sufficiently severe such that the failed spacecraft 
is no longer useful. In this case, it should be removed from the formation by deorbiting or by 
manual retrieval, If the failed spacecraft is replaced by a backup spacecraft, then it is necessary to 
reconfigure the formation, The control laws for steering the remaining active spacecraft from the 
old to the new formation requires separate consideration. This aspect will be discussed elsewhere. 

We note also that in the derivation of foregoing control laws, no constraints have been imposed 
on the magnitude of the control variables. I n the presence of bounded controls, one expects that 
the rate of decay of ||([ J Ej] i ,[^i]i)(f)|| and ||(5w,:,(5g. / )(f)|| to zero would be reduced by when one or 
more of the control variables takes on its extreme values. 

Finally, for a real mission, it is necessary to consider discrete-time versions of the proposed control 
laws. I n view of the limited fuel on-board, it is generally undesirable to have continuously acting 
controls. Therefore, the system response corresponding to the cent rol laws derived here serves as 
a basis for comparison between the idealized and the actual responses. 

6 Fleet Coordination 

For a fleet of spacecraft, one may require complete autonomy in each spacecraft in the sense that 
all the decisions for determining its future behavior are made on-board without the assistance 
of external agents. Although this approach provides enhanced operational reliability, it may not 
be cost effective since each spacecraft must contain all the essential hardware and software for 
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coordi nation and control. An alternative approach is to require each spacecraft to have only the 
basic hardware and software for attitude control and orbital nameuvering. The more complex tasks 
in fleet coordination and control are shared by all the spacecraft in the fleet. Moreover, some of 
the spacecraft may be equipped with special hardware and software to perform particular tasks for 
the entire fleet. 

The fleet coordination is achieved with the aid of an inter-spacecraft communication network (eg. 
radio or optical links). This network has the following basic functions: 

(i) Communicating the necessary data for fleet formation-keeping and relative attitude control; 

(ii) Linking the computers in the spacecraft to form a distributed computing network thereby 
increasing the computational capability of the fleet for more computational intensive tasks such as 
on-board interferometer data processing, 

In the realization of the first function, each fleet leader broadcasts its position and attitude 
with respect to a specified inertial frame, and the follower spacecraft broadcast their positions and 
velocities relative to their leader to achieve formation alignment. 


7 Simulation Studies 

Extensive simulation studies have been made to determine the performance of the proposed control 
laws for formation keeping and attitude regulation in the presence of actuator saturation, variations 
in spacecraft parameters, and loss of communication between spacecraft. Only typical results will 
be presented here. 

We assume that the leader of the spacecraft triad moves along an inclined circular orbit 0\ about 
the Earth with inclination angle (tt/2 - ipi nc ) and ascending node along the Y-axis. For convenience, 
weintroducea geocentric fixed cartesian coordinate frame T 0 with origin O at the Eart h's center 
along with a spherical coordinate system (r, 9, <f>) wit h orthonormal basis {e r , eg , }. Its motion 

in spherical coordinates is given by 

n(i) = To, 8i(t) = cos~ 1 {cos(<p inc )cos(8 0 - Uot)}, 

<j>i(t) = tan~ l {—tan(6 0 - uj 0 t)/sin(ip inc )}, (22) 

wherer Q is a given orbital radius, and u a =vW r o- Here, for simplicity, we have set the desired 
orbit al radius for all spacecraft to r Q . 

The desired motions for the second and third spacecraft correspond to two circular orbits with 
the same inclination angle (/x/2 - <pi nc = 8.2/i/180rad., but with ascending nodes at (r, Q,<f>) = 
(r Q , n/2, A(j)) and (r 0> (jl/2, - A <£>) respectively, where A <j) =9.5yu/180rad. We adopt the simplified 
control law ( 11" ) for formation keeping, where the deviation vector ^(f) is given by 

h 2 (t) = r 0 {(— cos(A<j>)sin (vw)cos(0(f)) + sin(A<f>)sin(6(t) )ex 

+ ( sin(A<j>)sin(<pi nc ) cos(9(t,)) + cos(A9)sin(9(t)) - sin(9 0 - w 0 f))ey 
+ (cos(ipi nc )cos(6(t) - A 9) - cos(0 o - 0J o t))ez- (23) 

where 9{t) = 19 0 - u Q t- A#. The deviation vector h$(t) has the same form as (21) except with A <j> 
replaced by - A <j>. Evidently, hi (t) satisfies d 2 hj, ( t)/dt 2 +u> 2 hi(t) = O for all t and % =2,3. To 
specify the desired attitude of thei-th follower, we introduce the (1-2-3) Euler angles (0,,^^^) 
corresponding to a rotation of © 7 ; about the X-axis followed by a rotation of ^ about the rotated 
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Y-axis, and a rotation of about the rotated Z-axis. The desired Euler angles for the followers 
are given by: 


02 (*) = © 3 (*) = m, *i(t) = *i(t) = Vine, 

*i(t) = A <f>, $g(i) = - A <t>. (24) 

Thus, the desired attitude of the i-th spacecraft can be expressed in terms of the following quater- 
nions: 


<724 (*) = {l + cos(A <f>)cos(<p inc ) - sin(Aip)sin(ip inc ) sin(d(t)) 

+ cos (0(f) ) (cos (A^) + cos(v?,; nc ))}2 /2; (25) 

<721 (f) = {sm(A0)sm(</?.; nc )cos(0(f)) + sin(8(t))(cos( A<j>) 

+ COS ((pine) )}/ ( 4<?24 (*) ) ; ( 26 ; 

922(*)-{a*n(¥>inc)(l + cos(A(j>)cos(6(i))) _ sin(A^)iin (fl(t))}/(4^(t)) ; ( 27 ) 

<723 (^) = {cos(A</>) sin(ipi nc )sin(8(t )) 

+ sin{A(j))cos(d(t )) + cos(v5inc))}/(4^ 4 (t)). (28) 

The quaternions corresponding to the desired attitude of the third follower have the same form as 
(23) except with A cj> replaced by - A+. We require every spacecraft to spin about its z-axis with 
constant angular speed ui s . Thus, the desired angular velocity for the i-th follower is given by 

ojf - u 0 cos((p inc )e x + L 0 o sin(<p inc )e z + u s e d iz = u> s ef x + u s e d z , (29) 

where {ef xl responds to the basis of the body coordinate system Tf associated with the 
i-th follower with the desired attitude. 

Figure 4 shows a typical time-domain response of the MS fleet with the simplified formation- 
keeping control law(ll” ) and attitude control law (19) in the presence of actuator saturation, The 
spacecraft parameter values used in the simulation study are given in Table 1. The corresponding 
time-domain response of the MS fleet with the [u c i]i term in (11") set to zero (to simulate the 
loss of communication between the MS) was also determined. The results do no differ significantly 
from those shown in Fig. 4. Next, the effect of inertia pertrubations on the time-domain response 
of the spacecraft triad was studied. It was found that the qualitative behavior of the response is 
essentially identical to that of the unperturbed case. 


8 Concluding Remarks 

In this paper, control laws for a spacecraft triad moving in formation have been derived using a 
simplified model for a rigid spacecraft. These control laws require the knowledge of the relative 
displacements and attitudes of the spacecraft and its neighbors. Simulation results based on a 
generic spacecraft model showed that the derived control law are effective in formation and rela- 
tive attitude alignment provided that the magnitude of the initial deviation from the desired state 
is sufficiently small so that collisions between the spacecraft do not occur. Finally, in this work, 
important factors such as data processing time-delay and time discretization arising in physical 
implementation have not been taken into consideration. Nevertheless, the results reveal the basic 
structure of the control laws and the required inter-spacecraft data required for their implemen- 
tation. Finally, the problems associated with the physical implementation of the control laws in 
terms of the state-of-the-art hardware and fuel consumption for control are not considered here, 
and they require further study. 
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Mi 

mass of M S 

10%. 

■fix 

moment of inertia about x-axis 

0.3646 kgm 2 . 

liy 

moment of inertia about y-axis 

0.2734%m 2 . 

Iiz 

moment of inertia about z-axis 

0.3125 kgm 2 . 

To 

desired orbital radius of MS 

7.13814 xl0 6 m. 

'o = \/aAo 

orbital angular speed of fleet leader 

0.001 rad. /sec. 

UJ S 

desired spin speed about z-axis 

0. 01 rad,./ sec. 

(fine — 7r/2 

inclination angle of reference orbits 

8.27 r/ 180rad. 

A<f 

azimuthal angle associated with the ascending node of reference orbits 0.2 rad.. 

A9 

MS separation angle 

7r/120 rad. 


Table 1: Values of microspacecraft and orbital parameters for simulation study, 
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List of Figures 


Fig. 1 Sketch of inertial and moving coordinate systems, 

Fig. 2 Exaggerated sketch of the reference orbits of spacecraft triad moving in formation along 
LEO about the Earth. 

Fig. 3 Time-domain response of the spacecraft triad with simplified formation-keeping control law 
(11") and attitude control law (19) with Ku = 0,5, = 2.0, K ui = 1.5, K qi = 

0.4, * = 2, 3, in the presence of actuator saturation, and with initial states: 

[£ 2 (0),J3 2 (0)] 2 = [5 2 -5, 00 O]; [w 2 (0)] 2 = [0.02 0,02 0.02]; 
g 2 (0) = [0.3 0.1 0.2 0.9272618]; 

[E 3 (0), £ 3 (0)] 3 = [-1 1 -1, 00 O]; [u*(0)] 3 =[-0.01 0.015 - 0.01]; 

<? 3 (0) = [-0.2 0.2 0.3 0.910433]; 

follower spacecraft 2 (solid curves); follower spacecraft 3 (dashed curves); Saturation levels: 
WfdjW < 1 N) || T ci j < 0.05 N.m., i = 2,3; j = x,y, z. 

Fig. 3a: Positional tracking errors (m) vs.time. 

Fig. 3b: Angular velocities (m/see) vs. time. 

Fig. 3c: Quaternions vs. time. 

Fig. 3d: Control forces (N) vs. time. 

Fig. 3e: Control torques (N.m) vs. time. 
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Figure 1. Sketch of inertial and moving coordinate systems. 
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LEADER 
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Figure 2. Exaggerated sketch of the reference orbits of 
spacecraft triad moving in formation along LEO about the 
Earth. 
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Figure 3. Time-domain response of the spacecraft triad with 
simplified formation-keeping control law (11”) and attitude 
control law (19) with K ll =0.5 > K 2I =3.0 1 K ai =1.5, ^=0.4, 
1=2,3, in the presence of actuator saturation, and with initial 
states: 

[E 2 (O),Ej(0)] 2 ={52 -5, 00 O]; 

(CO,(o)J, = (0.02 0.02 0.02); 

q 2 (0)= [0.3 0.10.20.927261 8]; 

[E 3 (0),E 3 (0)L = [-1 2-1, 0 0 0]; 

[oo 3 (0)n-0.01 0.015- 0.01]; 
q 3 (0) =[-0.2 0.20.3 0.910433]; 

The Follower spacecraft 2 is shown with solid lines and the 
follower spacecraft 3 is shown with dashed lines; Saturation 
levels: 

lf cij l < 1 N; lT cij i < 0.05 N.m, 1=2,3; j=x,y,z. 
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Figure 3a. Positional tracking errors (m) VS. time. 
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